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Abstract 



A search for pair production of second-generation scalar leptoquarks in the final state 
with two muons and two jets is performed using proton-proton collision data at y/s = 
7 TeV collected by the CMS detector at the LHC The data sample used corresponds 
to an integrated luminosity of 34 pb _1 . The number of observed events is in good 
agreement with the predictions from the standard model processes. An upper limit 
is set on the second-generation leptoquark cross section times fi 2 as a function of the 
leptoquark mass, and leptoquarks with masses below 394 GeV are excluded at a 95% 
confidence level for (6 = 1, where f> is the leptoquark branching fraction into a muon 
and a quark. These limits are the most stringent to date. 
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Several extensions of the standard model HH3 predict the existence of leptoquarks (LQ), hy- 
pothetical particles that carry both lepton and baryon numbers and couple to both leptons 
and quarks. Leptoquarks are fractionally charged and can be either scalar or vector particles. 
In order to satisfy constraints from flavour-changing neutral currents and rare pion and kaon 
decays ©13/ leptoquarks are restricted to couple to a single lepton-quark generation. 

In proton-proton collisions at the CERN Large Hadron Collider (LHC) the dominant mech- 
anisms for pair production of scalar leptoquarks are gluon-gluon fusion and ^-annihilation. 
The cross section depends on the strong coupling constant and the LQ mass and has been cal- 
culated at Next-to-Leading-Order (NLO) |8); the dependence on the Yukawa coupling A is 
negligible [8]. Leptoquarks decay to a quark and a charged lepton of the same generation with 
unknown branching fraction j5 and to a quark and a neutrino with branching fraction (1 — jS). 
In this analysis, we consider the decay of a second-generation leptoquark to a muon and a 
quark. 

Several experiments have searched for leptoquarks, but so far no evidence has been observed. 
A review of LQ phenomenology and searches can be found in (9l[l0|. The most recent limits 
from the D0 experiment at the Fermilab Tevatron collider exclude second-generation scalar 
leptoquarks with masses below 316 GeV for /3 = 1, based on proton-antiproton collisions at 
</s = 1.96 TeV El • 

This Letter describes a search for pair production of second-generation scalar leptoquarks with 
the CMS experiment using LHC proton-proton collisions at i/s = 7 TeV. The data sample used 
corresponds to an integrated luminosity of 34.0 ± 3.7pb _1 . 

The CMS detector, described in detail elsewhere |l2fl, uses a cylindrical coordinate system with 
the z axis along the counterclockwise beam direction The angles 9 and (p are the polar and 
azimuthal angles, respectively. Pseudorapidity is defined as rj = —In [tan (0/2)], where 9 is 
measured with respect to the +z-axis. The central feature of the CMS apparatus is a supercon- 
ducting solenoid, of 6 m internal diameter, providing a field of 3.8 T. Within the field volume 
are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter (ECAL) and the 
brass/ scintillator hadron calorimeter (HCAL). Muons are measured in gas-ionization detectors 
embedded in the steel return yoke. In addition to the barrel and endcap detectors, CMS has 
extensive forward calorimetry. The inner tracking system consists of a silicon pixel and strip 
tracker, providing the required granularity and precision for the reconstruction of vertices of 
charged particles having pseudorapidities \tj\ < 2.5. The ECAL and HCAL are used to measure 
the energies of photons, electrons, and hadrons within a region of \f]\ < 3.0. The three muon 
systems surrounding the solenoid cover a region \rj\ < 2.4 and are composed of drift tubes 
in the barrel region (\r/\ < 1.2), of cathode strip chambers in the endcaps (0.9 < \rj\ < 2.4), 
and of resistive plate chambers in both the barrel region and the endcaps (\tj\ < 1.6). Events 
are recorded based on a first-level trigger decision coming from either the calorimeter or muon 
systems. The final trigger decision is based on the information from all subsystems, which is 
passed on to the High Level Trigger (HLT), consisting of a farm of computers running a ver- 
sion of the reconstruction software optimized for fast processing. The signature of the decay 
of pair-produced second-generation leptoquarks studied here consists of two muons and two 
jets with high transverse momentum (pr)- Events are selected by a single muon trigger with- 
out isolation requirements and with lower pj thresholds dependent upon the instantaneous 
luminosity. The HLT efficiency for events with two high pj muons is always greater than 99%. 

The Monte Carlo (MC) signal events are generated using the PYTHIA [T3l| generator (version 
6.422) in the LQ mass range 250-600 GeV. The CTEQ6.6 [14J Parton Distribution Function (PDF) 
sets are used. The main background processes that can mimic the signature of the LQ signal 
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are Z/7* + jets, tt, VV (WW,ZZ,WZ), W + jets, and multijet events. The tt, VV, and muon- 
enriched multi-jets events are generated with MadGraph |fl5lH6} ; Z/7* + jets and W + jets 
events are generated with ALPGEN. In MADGRAPH and ALPGEN samples, parton showering 
and hadronization is performed with PYTHIA. 

Muons are reconstructed as tracks in the muon system that are matched to the tracks recon- 
structed in the inner tracking system. Muons are required to have pj > 30 GeV, \t]\ < 2.4. 
The muon relative isolation paramter is defined as the scalar sum of the pj of all tracks in 
the tracker and the transverse energies of hits in the ECAL and HCAL in a cone of AR = 
\J [At]) 2 + (A(p) 2 = 0.3 around the muon track, excluding the contribution from the muon it- 
self, divided by the muon pj. Muons are required to have a relative isolation parameter less 
than 0.05. Arj and A(p are the pseudorapidity and azimuthal angle differences between the 
muon track and other reconstructed tracks or hits in the calorimeter. To have a precise mea- 
surement of the transverse impact parameter of the muon relative to the beam-spot position, 
only muons with tracks containing more than 10 hits in the silicon tracker are considered. To 
reject muons from cosmic rays, the transverse impact parameter is required to be less than 
2 mm. In addition, the two muon candidates are required to be separated from each other by 
at least AR = 0.3 and at least one muon must be in the pseudorapidity region \rj\ < 2.1. The 
efficiency of selecting dimuon events is 61-70% for the LQ mass range of 200-500 GeV. 

Jets are reconstructed using the anti-kj |[T7| algorithm with a distance parameter R = 0.5 and 
are required to have pj > 30 GeV and \rj\ < 3.0. Jet-energy-scale corrections derived from 
MC simulated events are applied to establish a relative uniform response in \] and an absolute 
uniform response in pj. A residual jet energy correction is derived from data by looking at the 
balance in pj in dijet events, and it is applied to jets in data. 

Additional selection requirements are placed on two variables, which are effective at discrimi- 
nating the LQ signal from the major sources of background. The first is the dimuon invariant 
mass, Mpp. The second variable, St, is defined as the sum of the magnitudes of the pr of the 
two highest pr muons and the two highest pj jets. The two muons in the signal events come 
from the decays of two high-mass particles, and they tend to form a large invariant mass. Thus, 
events are selected if M w > 115 GeV. This helps to reduce the contribution from Z/7* + jets 
processes, which is one of the largest backgrounds. In addition, the LQ pair is expected to have 
a large St- The lower threshold on St is optimized for different LQ mass hypotheses by using a 
Bayesian approach |Tl8llT9l to minimize the expected upper limit on the LQ cross section in the 
absence of an observed signal. The St cut helps to further reduce background sources, most 
noticeably tt. The optimal St threshold values for each mass hypothesis are given in Table [l] 
While the LQ signal is expected to peak in the mass distribution of the p-jet pairs, we find that 
the St variable gives sufficient power of discrimination in the range of LQ masses considered. 
The p-jet mass distribution would nevertheless be important to establish the signal in case an 
excess is observed. 

The contribution from tf is estimated with the MC sample, using uncertainties determined from 
data [20 [. The contribution from W + jets is negligible once the full event selection is applied. 
The small contribution from VV is estimated from MC. The multijet background is found to 
be negligible using a control data sample of same-sign dimuon events. The background from 
Z/ 7* + jets is determined by comparing Z/ 7* + jets events from data and MC samples in two 
different regions: at the Z boson peak, 80 < < 100 GeV, and in the high-mass region, 

M }l pi > 115 GeV. In the low-mass region, the ratio of data to MC events (Rl) is determined to be 
Ri = 1.28 ± 0.14 after selecting two muons and two jets with px > 30 GeV, and a preliminary 
requirement of St >250 GeV. This rescaling factor is applied to the number of Z/ 7* + jets MC 
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events in the high-mass region after the full selection. 

Reasonable agreement between data and MC predictions is observed at all selection levels. The 
dimuon invariant mass is shown in Fig. [T] (left) after the initial selection of muons and jets with 
pr > 30 GeV and a preliminary requirement of St > 250 GeV. The M w distribution in data is 
consistent with the expected SM background prediction. The St distribution is also shown in 
Fig. [I] (right) after the initial selection of muons and jets with pj > 30 GeV and the additional 
requirement of M flfl > 115 GeV. 
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Figure 1: The distribution of M™ (left) after requiring at least two muons and at least two 
jets with pt > 30 GeV and St > 250 GeV arid the distribution of St (right) after requiring at 
least two muons and at least two jets with pt > 30 GeV and M w > 115 GeV. The Z/ 7* — > 
+ jets contribution is rescaled by the factor Rl = 1.28. Other backgrounds correspond to 
VV, W + jets, and multijet processes. Uncertainties are statistical. 



The event yields from data, expected LQ signal (for several mass hypotheses), signal selection 
efficiency times acceptance, and expected standard model backgrounds are summarized in 
Tabled] 

Table 1: The data event yields in 34.0 pb -1 for different leptoquark mass hypotheses, together 
with the optimized St threshold values (in GeV) for each mass, background predictions, num- 
ber of expected LQ signal events (S), and signal selection efficiency times acceptance (€$). Mlq 
and St values are listed in GeV. The Z/ 7* — > fiji + jets contribution is rescaled by the factor 
Rl = 1.28. Other backgrounds correspond to VV, W + jets, and multijet processes. Uncertain- 
ties are statistical. 



M LQ 


MC Signal Samples 


Monte Carlo Background Samples 


Events 


Obs./Exp. 


(S T Cut) 


Selected 


Acceptance 




Selected Events in 




in 


95% C.L. 


[GeV] 


Events 


x Efficiency 


tt + jets 


Z/7* + jets 


Others 


All 


Data 


u.l. on a [pb] 


200 (S T > 310) 


160±20 


0.388±0.003 


3.9±0.1 


4.08±0.07 


O.lrtO.Ol 


8.1±0.1 


5 


0.661/0.444 


225 (S T > 350) 


89±9 


0.421±0.003 


2.7±0.1 


2.99±0.05 


0.07±0.01 


5.8±0.1 


3 


0.525/0.342 


250 (S T > 400) 


51±5 


0.437±0.003 


1.6±0.08 


1.92±0.04 


0.051±0.009 


3.57±0.09 


3 


0.422/0.370 


280 (S T > 440) 


28±3 


0.467±0.003 


0.98±0.06 


1.53±0.03 


0.038±0.008 


2.55±0.07 


3 


0.352/0.375 


300 (S T > 440) 


21±2 


0.518±0.004 


0.98±0.06 


1.53±0.03 


0.038±0.008 


2.55±0.07 


3 


0.318/0.339 


320 (S T > 490) 


14±1 


0.509±0.004 


0.55±0.04 


1.12±0.02 


0.019±0.005 


1.69±0.05 


2 


0.287/0.302 


340 (S T > 530) 


9±1 


0.508±0.003 


0.34±0.04 


0.79±0.01 


0.01±0.004 


1.14±0.04 


1 


0.261/0.246 


400 (S T > 560) 


4.0±0.4 


0.578±0.004 


0.27±0.03 


0.67±0.01 


0.01±0.004 


0.95±0.03 





0.220/0.219 


450 (S T > 620) 


1.9±0.2 


0.600±0.004 


0.16±0.03 


0.49±0.01 


0.006±0.003 


0.66±0.03 





0.198/0.153 


500 (S T > 700) 


0.9±0.1 


0.602±0.004 


0.08±0.02 


0.277±0.006 


0.003±0.002 


0.36±0.02 





0.179/0.152 



Several sources of systematic uncertainties are considered in this analysis. The uncertainty on 
the integrated luminosity is taken as 11% I|2T1 . A 5% systematic uncertainty is assigned to the 
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jet-energy-scale (JES) [22 j of each jet. A smaller, ~ 1% systematic uncertainty comes from the 
muon momentum scale. The 300 GeV LQ signal efficiency changes by 2% and 1% due to JES 
and muon momentum scale uncertainties, respectively. The effect of the muon momentum 
scale uncertainty on the total background is estimated to be < 0.5%. The JES contributes 2% to 
the estimate of the Z/ 7* + jets background described above and 15% to the estimate of the VV 
background from MC. The statistical uncertainty on the value of Ri after a preselection require- 
ment (St > 250 GeV), 11%, is used as an uncertainty on the estimated Z/ 7* + jets background. 
A 41% systematic uncertainty is taken from the CMS measurement of the tt production cross 
section pOfl and assigned to the estimate of the tt background; it includes the effect of JES on 
the estimate of the tt background. The effect of jet energy and muon momentum resolution on 
expected signal and backgrounds is found to be negligible. A 5% systematic uncertainty per 
muon is assigned due to differences in reconstruction, identification, trigger, and isolation effi- 
ciencies between data and MC 1123 1 , resulting in a 10% uncertainty on the efficiency of selecting 
events with two muons both for the signal and background processes. A theoretical uncertainty 
on the LQ signal production cross sections due to the choice of renormalization/ factorization 
scales has been calculated by varying the scales between half and twice the LQ mass, and is 
found to 14-15% for LQ masses between 200 and 500 GeV. The 90% C.L. PDF uncertainties on 
LQ cross section have been obtained using the CTEQ6.6 |14| PDF error set following a standard 
prescription and have been found to vary from 8 to 22 % for leptoquarks in the mass range of 
200 - 500 GeV (8]|. The effect of PDF uncertainties is less than 0.5% on signal acceptance. The 
PDF uncertainties are not considered for background sources with uncertainties determined 
from data. The systematic uncertainties, their magnitude, and the relative impact on the num- 
ber of signal and background events are summarized in Table |2j 

Table 2: Systematic uncertainties and their effects on number of signal and background events 



Systematic 


Magnitude 


Effect on 


Effect on 


Uncertainty 




Signal 


Background 


JES 


5% 


2% 




JES & Data Backgr. Est. 






25% 


Muon Momentum Scale 


1% 


1% 


< 0.5% 


Muon Pair Reco/ID/Iso 


10% 


10% 


< 0.05% 


Integrated Luminosity 


11% 


11% 




Total 




15% 


25% 



One candidate event survives the full selection criteria corresponding to a leptoquark mass 
hypothesis of 340 GeV, and no candidates survive for criteria corresponding to masses greater 
than 450 GeV. An upper limit on the LQ cross section is set using a Bayesian method |Tl8l H"9l 
with a flat signal prior. A log-normal probability density function is used to integrate over the 
systematic uncertainties. Using Poisson statistics, a 95% confidence level (C.L.) upper limit is 
obtained on a x f} 2 . This is shown in Fig. [2] together with the NLO predictions for the scalar 
LQ pair production cross section. The 95% C.L. exclusion on /6 as a function of LQ mass is also 
shown in Fig. [2] The systematic uncertainties reported in Tableware included in the calculation 
as nuisance parameters. With the assumption that /3 = 1, second-generation scalar leptoquarks 
with masses less than 394 GeV are excluded at 95% C.L., 78 GeV higher than the limit set 
at the D0 Experiment at the Tevatron flTTfl. This is in agreement with the expected limit of 
394 GeV. The corresponding observed limit on cross section is 0.223 pb. If the lower edge of the 
theoretical a x /S 2 curve is used, the observed (expected) limit on LQ mass is 380 (378) GeV and 
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the observed limit on cross section is 0.225 pb. 
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Figure 2: (Left) The expected and observed 95% C.L. upper limit on the scalar leptoquark pair 
production cross section multiplied by /3 2 as a function of the LQ mass, together with the NLO 
theoretical cross section curve. The shaded band on the theoretical values includes CTEQ6.6 
PDF uncertainties and the error on the leptoquark production cross section due to renormaliza- 
tion and factorization scale variation by a factor of two. The shaded region is excluded by the 
current D0 limits fU\. (Right) The minimum /3 for 95% C.L. exclusion of the leptoquark hy- 
pothesis as a function of leptoquark mass. The observed limit and corresponding uncertainty 
band is obtained by considering the observed upper limit and theoretical branching ratio and 
its uncertainty in the left-hand figure. Note: The shaded area excluded by the D0 experiment 
was determined with combined information from the decay channel with two muons and two 
jets and the decay channel with one muon, missing transverse energy and two jets. 



In summary a search for pair production of second-generation scalar leptoquarks decaying to 
two muons and two jets has been performed using 7 TeV pp collision data corresponding to 
an integrated luminosity of 34.0 pb -1 . The number of observed candidate events agrees well 
with the number of expected standard model background events. A Bayesian approach that 
includes the treatment of systematic uncertainties as nuisance parameters is used to set limits 
on the LQ cross section times /5 2 as a function of LQ mass. At 95% C.L., the pair production of 
second-generation scalar leptoquarks with masses below 394 GeV is excluded for j6 = 1, where 
j6 is the leptoquark branching fraction into a muon and a quark. This is the most stringent limit 
to date on the existence of second-generation scalar leptoquarks. 
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